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SUMMARY 


Supersonic  profiles  of  minimum  pressure  drag  for  a  given  thickness 
ratio  and  for  a  given  area  have  been  determined  with  the  use  of  a  non¬ 
linear  pressure  relation  and  are  compared  with  minimum- drag  profiles 
found  by  linearized  theory.  The  results  show  that  the  profiles  are 
determined  with  sufficient  accuracy  by  linear  theory  over  the  entire 
supersonic  Mach  number  range  since  the  drag  coefficients  for  these 
profiles  are  only  slightly  higher  than  those  for  optimum  profiles  deter¬ 
mined  by  nonlinear  theory.  Linear  theory  appears  to  be  adequate  for 
determining  profiles  of  minimum  drag  for  other  auxiliary  structural 
conditions  since  moderate  deviations  from  the  optimum  shape  have  only 
a  small  influence  on  the  pressure  drag. 

The  parameters  determining  the  airfoil  shape  for  a  given  thickness 
ratio  found  by  both  the  linear  and  nonlinear  theory  are  presented  in 
graphs  as  a  function  of  the  base  pressure  coefficient.  With  the  use  of 
these  results,  the  optimum  profiles  for  any  stream  Mach  number  and  thick¬ 
ness  coefficient  are  readily  determined.  A  comparison  of  the  pressure 
drag  coefficients  for  optimum  profiles  determined  by  linear  and  nonlinear 
theory  is  presented  for  the  Mach  number  range  from  1.5  to  10.0.  In 
addition,  several  optimum  profiles  for  a  given  area  have  been  calculated 
by  both  the  linear  and  nonlinear  theory. 


INTRODUCTION 


Drougge  (reference  l)  has  determined  the  airfoil  section  shape  for 
minimum  pressure  drag  at  supersonic  speeds  subject  to  such  auxiliary 
conditions  as  given  bending  and  torsional  stiffness.  These  calculations 
were  made  by  using  the  linearized  expression  for  the  pressure  coefficient; 
the  effect  of  a  base  was  not  considered.  Recently,  Chapman  (reference  2) 
has  shown  that  the  section  shape  for  minimum  pressure  drag  as  determined 
by  linearized  theory  may  have  a  blunt  trailing  edge.  The  use  of  linear¬ 
ized  theory  for  determining  optimum  profiles  facilitates  the  mathematical 
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development;  however,  the  results  are  subject  to  question  particularly 
at  high- supersonic  Mach  numbers. 

w 

The  purpose  of  the  present  paper  is  to  compare  the  section  shapes 
for  minimum  pressure  drag  (subject  to  certain  auxiliary  conditions) 
determined  by  linear  and  nonlinear  theory  In  order  to  estimate  the 
errors  introduced  by  the  linearized  form  of  the  pressure  coefficient 
and  to  determine  its  range  of  validity  for  calculations  of  this  nature. 

For  this  purpose,  it  was  considered  sufficient  to  examine  two  problems. 

The  problems  chosen  were  the  determination  of  the  profile  for  minimum 
drag  for  a  given  thickness  ratio  and  the  determination  of  the  profile 
for  minimum  drag  for  a  given  area. 

The  nonlinear  form  of  the  expression  for  the  pressure  coefficient 
used  in  the  present  analysis  is  derived  in  reference  3  where  it  is 
shown  to  be  in  excellent  agreement  with  the  exact  expression  for  stream 
Mach  numbers  greater  than  1.5-  The  variation  of  base  pressure  coeffi¬ 
cient  with  stream  Mach  number  was  assumed  in  order  to  determine  actual 
airfoil  shapes.  This  base -pres sure  curve  was  based  in  part  upon  experi¬ 
mental  data  and  the  known  variation  of  vacuum  pressure  coefficient  with 
Mach  number. 
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NOTATION 
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00 


chordwise  distance  from  leading  edge 
airfoil  ordinate 

4 

airfoil  chord;  also,  — - ■ 

ratio  of  specific  heat  at  constant  pressure  to  specific  heat 
at  constant  volume 


stream  Mach  number 


P 


pressure  coefficient 
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3 


t  thickness  coefficient 


3  flow  deflection  angle 

0  =  m3 

/ 

c^  drag  coefficient 

A  airfoil  area 

Subscripts: 

front  surface 
rear  surface 

base  \ 

vacuum 

refers  to  double -wedge  profile 
optimum 
trailing  edge 

ANALYSIS  AND  DISCUSSION 
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The  pressure  drag  of  an  airfoil  is  the  drag  due  to  the  normal' forces 
acting  over  the  airfoil  surface.  For  a  profile  with  a  blunt  trailing 
edge  the  base  drag  must  be  added  to  the  drag  due  to  the  normal  forces 
over  the  forward  surfaces  to  give  the  total  pressure  drag.  With  the 
surface  pressure  coefficient  denoted  by  P^  the  airfoil  shape  by  •  Y(X), 
the  chord  by  c,  the  base  height  by  Y^,  and  the  base  pressure  coef¬ 
ficient  by  P^,  the  pressure  drag  coefficient  c^  of  a  thin  symmetric 
two-dimensional  airfoil  at  zero  angle  of  attack  is  given  approximately  by 
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X 


or  with  x  =  — ,  y  -  — ,  and  =  — ,  the  pressure  drag  may  he  written  as 
C  C  c 


c,  =  2  /  P  ^  dx 

d  |  dx 


Vb 


(1) 


Profile  for  minimum  drag  for  a  given  thickness  ratio.-  In  order  to 
determine  the  profile  for  minimum  drag  for  a  given  thickness  ratio,  it 
is  required  to  determine  the  function  y(x)  and  the  "base  height  y-^ 
which  make  c^  a  minimum.  For  supersonic  flow,  to  a  high  order  of 

approximation,  the  pressure  coefficient  is  a  function  only  of  the  local 
dv 

slope  ^  and  the  stream  Mach  number.  By  use  of  the  calculus  of  varia¬ 


tions,  any  integral  of  the  form 


F(-^L]dx  can  he  shown  to  have 
\dx/ 


a  stationary  value  for  =  Constant  (reference  4,  for  example).  Thus, 


since  the  integrand  of  equation  (l)  is  a  function  only  of  the  slope 
the  profile  of  minimum  drag  has  plane  surfaces.  Then,  with  $  =  ^  and 
the  notation  of  figure  1,  the  drag  coefficient  may  he  written  as 


c  =  2x1*1P1  +  2(1  -  xJflgPj,  -  ybPb 


(2) 


where  x^_  is  the  location  of  maximum  thickness  and 

P1  »  P(*l)  P2  »  P(*2) 


Then  with 


Xb 

t 


where  t  is  the  thickness  coefficient,  equation  (2)  may  be  written  as 


hP2  +  (h  -  l)Pb 


I 

f 


(3) 


naca'tn  2623 


5 


Reference  3  shows  that 


where  Mm  is  the  stream  Mach  number  and  7  is  the  ratio  of  specific 
heats  at  constant  pressure  and  constant  volume,  is  in  excellent  agree¬ 
ment  with  the  pressure  coefficients  given  by  the  exact  shock  or  expan¬ 
sion  relations  for  1.5  =  M^,  <  <».  It  may  be  noted  that  this  expression 
approaches  both  the  linearized  value  and  the  hypersonic  value  in  the  appro¬ 
priate  Mach  number  ranges.  With  6  =  md  as  the  independent  variable, 
equation  (4)  becomes 

P  =  K(0) 

'  '  m2c 

where 


Then,  with  these  relations  substituted  into  equation  (3),  the  drag 
coefficient  becomes 


cd  =  -#1  “  ^2  +  (*  “  D|  (5) 

m  e  L_  ^  S3 


where 


i 


e 


mt 


mt  h 
2  1  -  xx 
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Values  of  and  h  are  required  which  make  the  drag  a  minimum. 

Thus,  solutions  are  required  for  the  equations 


^ca 

dxj_ 


0 


for  and  h . 


The  base  pressure  coefficient  is  dependent  upon  Mach  number,  airfoil 
shape,  Reynolds  number,  the  nature  of  the  boundary  layer,  and  so  forth. 

Of  these  the  Mach  number  is  by  far  the  most  important  parameter.  Since 
the  base  pressure  depends  to  some  extent  on  the  base  height,  a  term 
C>P, 

involving  should  be  included  in  the  second  equation  used  to  deter¬ 


mine  the  airfoil  of  minimum  drag.  The  omission  of  this  term  however  is 
of  little  consequence.  Its  effect  is  to  modify  slightly  the  effective 
value  of  the  base  pressure.  For  turbulent  boundary  layers  at  high 
Reynolds  numbers  this  change  is  small  and  the  change  in  airfoil  shape 
due  to  the  neglect  of  this  term  is  of  little  consequence.  If  the  term 
Bpb 

is  omitted,  the  equations  defining  the  airfoil  of  minimum  drag  are 


iKi 


+  e22K2 


=  0 


(6) 


■5F  =  (“K2  +  e2K2*  +  I  m2pb)  =  0  (7) 

where 

K'  =  &  =  - _ - _ 

e2/e2  +  c2 

From  the  solutions  of  equations  (6)  and  (7)  for  0 ^  and  02  as 
functions  of  (m^  -  l)^,  the  location  of  maximum  thickness  and  the 
trailing-edge  thickness  are  determined  from 


NACA  TN  2623 


7 


•\ 


ei02 


“ 


The  lower  limit  for  nrt  insures  that  the  rear  upper  and  lower  surfaces 
can  intersect  only  at  the  trailing  edge.  The  upper  limit  for  mt 
insures  that  the  location  of  maximum  thickness  will  lie  between  the 
leading  and  trailing  edges  of  the  airfoil. 

It  is  important  to  note  that  0-^  and  ©2  are  functions  only  of 
(Moo2  -  ljPt).  Therefore,  x-|_  and  yb/t  are  functions  only  of  t^M^2  -  1 
and  (Mo.2  -  l^Pfc.  The  number  of  parameters  has  thus  been  reduced  from 
three  to  two.  From  equation  (8)  the  trailing-edge  thickness  is  zero 
^1^2 

for  mt  =  -2- - — ;  thus,  the  profile  of  minimum  drag  has  a  sharp 

dl  ' 

trailing  edge  for  this  condition.  For  profiles  with  a  sharp  trailing 
edge,  only  equation  (6)  need  he  solved  and  the  location  of  maximum 

thickness  then  is  a  function  only  of  t^M^2  -  1.  For  mt  =  £©-]_,  the 
the  position  of  maximum  thickness  for  minimum  drag  is  at  the  trailing 
edge.  For  values  of  mt  >  20-^ ,  no  mathematical  minimum  for  the  drag 
exists. 

The  solution  of  the  equations  defining  the  profile  of  minimum 
pressure  drag  is  especially  simple  for  linearized  flow  (reference  2, 
for  example).  The  approximate  expression  for  the  pressure  coefficient 
given  by  linearized  theory  is 


and  the  equations  which  determine  the  profile  of  minimum  drag  are 
8c 


8x- 


=  =  P  or  @i  =  -6: 


(9) 
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8cd 

"5h 


m2Pi 


=  0 


or  e2 


2-p 

m  Pb 

“T~ 


(10) 


Then,  the  location  of  maximum  thickness  and  the  trailing-edge  thickness 
are  found  to  be 


xi  = 

yb 

t 


2t 

mPv 


2t 


+  2  =  2- 


(■ 


_2-p  ' 

"  pb  <  <  "A 

— Jj—  -  Bit  -  -  2  , 


For  a  profile  with  a  sharp  trailing  edge 


(n) 


the 


familiar  result  that,  for  minimum  drag,  the  location  of  maximum  thickness 
is  at  the  midchord  is  found  as  the  solution  of  equation  (9)  for  y^  =  0. 


The  solution  of  equations  (6)  and  (7)  is  shown  in  figure  2  where 
and  62  are  given  as  functions  of  ^M^  -  l^P-^.  The  dashed  vertical 

line  at  -  l)Pfo  equal  to  1.^28  is  the  limiting  value  of  the 

vacuum  pressure  coefficient  for  infinite  Mach  number.  The  solutions  for 
linearized  flow  are  also  shown  in  figure  2  for  comparison.  It  may  be 
noted  that  the  values  of  for  linearized  flow  are  in  good  agreement 

with  the  results  of  this  analysis,  whereas  shows  appreciable  devia¬ 

tion.  Curves  of  t [/m^  -  1  are  also  shown  as  functions  of  (m^  - 

in  figure  2.  The  upper  set  of  curves  gives  the  value  of  -  1  for 

which  the  optimum  profile  has  the  maximum-thickness  location  at  the 

trailing  edge  and  the  lower  set  of  curves  gives  the  value  of  t{u^  -  1 
for  which  the  optimum  profile  has  a  sharp  trailing  edge.  At  a  given  value 

of  (m^2  ~  true  minimum-drag  profiles  do  not  exist  for  values  of 

t^Moo2  -  1  greater  than  the  values  given  by  the  upper  curve.  The  profile 
of  least  drag,  however,  under  those  conditions  is  one  which  has  the _ 

maximum-thickness  location  at  the  trailing  edge.  For  values  of 
less  than  those  of  the  lower  curve,  the  profile  of  minimum  drag  has  a 
sharp  trailing  edge  and  the  position  of  maximum  thickness  for  minimum  drag 
may  be  determined  from  figure  3*  Profiles  for  minimum  pressure  drag  are 
readily  determined  with  the  use  of  figures  2  and  3-  For  a  given  base 
pressure  coefficient  and  Mach  number,  values  of  6^  and  @2  are  found 
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Q  Q  _ 

from  figure  2.  Then  for  -2g-  -  1  =  29±,  the  position 

of  maximum  thickness  and  the  trailing-edge  thickness  are  found  from 
equations  (8)  for  a  given  thickness  coefficient.  For  values  of 


_  Q  Q 

t(/M  2  -  1  =  -2  1  ,  the  profile  has  a 

61  "  °2 

position  of  maximum  thickness  for  minimum 
from  figure  3 • 


sharp  trailing  edge  and  the 
drag  is  determined  directly 


In  order  to  show  more  clearly  the  variation  of  the  shape  of  the 
airfoil  for  minimum  drag  with  Mach  number,  a  number  of  profiles  have 
been  calculated  for  the  base  pressures  shown  in  figure  4.  The  base 
pressure  coefficients  in  this  figure  are,  in  part,  based  upon  some 
knowledge  of  the  base  pressures  for  turbulent  boundary  layers  and,  in 
part,  upon  the  variation  of  the  vacuum  pressure  coefficient  with  Mach 
number.  For  comparison  the  vacuum  pressure  coefficient  Pv  is  also 
shown.. 

In  figure  5  the  position  of  maximum  thickness  and  the  trailing- 
edge  thickness  for  minimum  drag  are  presented  as  functions  of  Mach 
number  for  the  base  pressure  coefficients  of  figure  4.  Here  it  may  be 
noted  that,  for  a  given  thickness  ratio,  the  minimum-drag  airfoil  has 
a  sharp  trailing  edge  for  the  lower  Mach  numbers  and  at  higher  Mach 
numbers  the  trailing  edge  is  blunt.  Further,  the  Mach  number  at  which 
the  minimum-drag  profile  first  has  a  blunt  trailing  edge  is  higher  for 
the  thinner  airfoils.  For  a  given  thickness  ratio,  the  position  of 
maximum  thickness  moves  rearward  with  increasing  Mach  number  until  it 
is  located  at  the  trailing  edge. 


Also  shown  in  figure  5  are  the  optimum  profiles  for  the  6-  and 
10-percent-thick  airfoils  calculated  from  the  linearized  equations  (ll). 
For  the  6-percent -thick  airfoil  the  position  of  maximum  thickness  remains 
fixed  at  the  midchord  up  to  a  stream  Mach  number  of  approximately  5 •  At 
Mach  numbers  greater  than  5,  the  position  of  maximum  thickness  moves 
rearward  and  the  trailing-edge  thickness  increases  in  a  manner  similar 
to  that  determined  from  the  nonlinear  equations.  For  the  10-percent- 
thick  profile  the  position  of  maximum  thickness  and  the  trailing-edge 
thickness  show  an  erratic  behavior  at  the  low  Mach  numbers,  that  is, 
at  low  values  of  the  base  pressure  coefficient.  At  a  stream  Mach  number 
of  1.5  the  optimum  position  of  maximum  thickness  is  at  the  0.56  chord. 

At  Mach  numbers  from  1.5  to  2.1  the  maximum-thickness  location  moves 
forward,  and  at  higher  Mach  numbers  it  moves  rearward  in  a  manner  similar 
to  that  determined  from  the  nonlinear  equations.  In  general,  the  linear¬ 
ized  theory  predicts  the  position  of  maximum  thickness  fairly  accurately 
over  a  wide  Mach  number  range  but  does  not  predict  the  optimum  trailing- 
edge  thickness  so  well. 
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In  order  to  illustrate  the  effect  of  the  airfoil  geometry  on  drag, 
the  drag  coefficients  of  the  optimum  profiles  have  been  calculated  for 
the  6-  and  10-percent -thick  airfoils.  In  figure  6  the  ratio  of  the 

optimum  drag  coefficient  c,  to  the  drag  coefficient  of  a  double- 

Qopt 

wedge  airfoil  c^;  is  shown  as  a  function  of  the  stream  Mach  number 

for  the  profile  determined  by  both  the  nonlinear  and  linear  theory.  For 
the  10-percent -thick  airfoil  the  drag  presented  is  not  a  minimum  above 
a  Mach  number  of  approximately  6.  The  drag  curve  above  this  Mach  number 
is  that  for  an  airfoil  with  the  maximum  thickness  located  at  the  trailing 
edge.  Here  it  may  be  noted  that  the  optimum  profiles  have  substantially 
lower  drag  than  the  double-wedge  profiles,  particularly  at  high  Mach 
numbers.  Further,  the  drag  reduction  is  considerably  greater  for  the 
10-percent-thick  airfoil  than  for  the  6-percent -thick  airfoil.  The 
dashed  curves  of  figure  6  correspond  to  the  drag  coefficients  for  the 
profiles  found  by  linear  theory.  In  calculating  the  drag,  however,  the 
nonlinear  expression  for  the  pressure  coefficient  was  used;  that  is, 
only  the  geometry  of  the  airfoils  was  calculated  by  linear  theory.  Here 
it  may  be  noted  that  the  drag  for  the  profiles  determined  by  linear 
equations  is  only  slightly  higher  than  the  drag  of  the  optimum  profiles 
found  with  the  nonlinear  relations.  The  largest  difference  in  drag 
between  the  6-percent -thick  profiles  determined  from  the  nonlinear  equa¬ 
tions  and  the  profiles  determined  from  the  linear  equations  is  less  than 
4  percent;  for  the  10-percent -thick  profiles  the  differences  in  drag  are 
even  less. 

Profile  for  minimum  drag  for  a  given  area. -  The  problem  of  determining 
the  profile  of  minimum  drag  and  satisfying  a  given  structural  condition 
is  an  isoperimetric  problem  of  the  calculus  of  variations.  The  equations 
for  determining  the  supersonic  profile  of  minimum  drag  for  a  given  struc¬ 
tural  condition  have  been  developed  in  reference  2.  The  equations  for 
determining  the  airfoil  of  minimum  drag  for  a  given  area  (or  torsional 
strength)  are 


J 


y(o)  =  0 
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where  y'  =  ^  X  is  a  constant  to  he  determined  from  the  solution, 

A  is  the  area,  and  the  subscript  TE  refers  to  the  values  at  the 
trailing  edge  of  the  airfoil.  The  solution  of  these  equations  for  the 
linearized  form  of  the  pressure  coefficient  gives  the  profile  shape  as 


-itK) 

where  is  the  position  of  maximum  thickness  and  is  given  by 


1  - 


xi = 


mP, 

b 

12A 


1  - 


mPb 

IT 


and  the  thickness  coefficient  is  given  by 


t 


The  solution  of  equations  (12)  for  the  nonlinear  form  of  the  pressure 
coefficient  was  obtained  by  an  iterative  procedure. 


Figure  7  presents  a  comparison  of  the  profiles  determined  by  linear 
and  nonlinear  theory  for  A  =  0.05  for  several  Mach  numbers  and  for 
the  base  pressure  coefficients  given  in  figure  h.  The  linearized  theory 
gives  a  location  of  maximum  thickness  farther  forward  and  a  smaller 
trailing-edge  thickness  and  thickness  ratio  than  given  by  nonlinear  theory. 
These  differences  in  geometry  of  the  profiles  determined  from  linear 
and  nonlinear  theory  follow  the  same  trend  as  for  the  airfoils  of  a  given 
thickness  ratio;  this  trend  may  be  expected  for  other  auxiliary  structural- 
conditions  as  well.  The  drag  of  the  profiles  found  by  linear  and  non¬ 
linear  theory  (equation  (4)  was  used  for  computing  the  pressures  for  the 
evaluation  of  the  drag  in  each  case)  differed  by  less  than  2  percent  for 
Mach  numbers  from  2  to  10. 


CONCLUDING  REMARKS 


The  results  show  that  the  profile  shape  for  minimum  drag  for  a  given 
thickness  ratio  or  for  a  given  area  is  determined  with  sufficient  accuracy 
by  linear  theory  over  the  entire  supersonic  Mach  number  range  since  the  drag 
coefficients  for  these  profiles  are  only  slightly  higher  than  those  for 
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optimum  profiles  determined  "by  nonlinear  theory.  It  would  appear  that 
linear  theory  should  also  be  adequate  for  determining  profiles  for 
minimum  drag  for  other  auxiliary  structural  conditions  since  moderate 
deviations  from  the  optimum  shape  have  only  a  small  influence  on  the 
pressure  drag  coefficient. 

It  appears  that,  when  airfoils  with  finite  trailing-edge  thicknesses 
are  considered,  the  linearized  theory  may  be  used  for  determining  profiles 
of  minimum  drag  (although  not  the  drag  itself)  at  least  up  to  Mach  numbers 
of  10. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  November  l4,  1951 
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functions  of  "base  pressure  and  the  limits  on  -  1  for  which 

the  solutions  are  valid. 
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Figure  4.-  Vacuum  pressure  coefficient  and  "base  pressure  coefficient  as 

functions  of  stream  Mach  number. 
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of  double-wedge  airfoil. 
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